Interfacial Energy and Fine Defect Structures 
for Incoherent Films 

Paolo Cermelli 

Dipartimento di Matematica, Universita di Torino, Italy 

Morton E. Gurtin 
Department of Mathematical Sciences, Carnegie-Mellon University, Pittsburgh, PA, U.S.A. 

Giovanni Leoni 
Dipartimento di Scienze e Tecnologie Avanzate, Universita del Piemonte Orientale, Italy 






o 

> 

O 
en 

o 
a^ 

-(— • 



o 



% 



This note summarizes recent results ||l| in which modern 
techniques of the calculus of variations are used to obtain 
qualitative features of film-substrate interfaces for a broad 
class of interfacial energies. In particular, we show that the 
existence of a critical thickness for incoherency and the for- 
mation of interfacial dislocations depend strongly on the con- 
vexity and smoothness of the interfacial energy function. 
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I. INTRODUCTION 

The structure of an interface separating crystalline 
solids is strongly affected by the competition between the 
relaxation of the bulk crystals to their respective equilib- 
rium configurations, and the tendency of the interface to 
mantain the exact matching (coherency) of the atoms of 
the two solids across the interface (cf. Leo & Hu [g|J|] 
and the review paper Q). When the stresses due to the 
deviation from equilibrium of the bulk material reach a 
suitable threshold, interfacial dislocations appear that re- 
lax the bulk stresses and an extreme situation may be 
reached in which all regularity of the atomic bonding at 
the interface is lost. 

We restrict attention to a two-dimensional framework 
with corresponding cartesian coordinates {x, y) and basis 
(i,j). We take the layer to be infinite in the x-direction, 
denote by h its height in the y-direction, and assume that 
2/ = represents the interface between the layer and the 
substrate, which we assume to be rigid. Letting u(a;, y) 
denote the displacement of the layer, e = ^(Vu -I- Vu^) 
the infinitesimal strain, and 

eo = eoi ® i, 

with eo > 0, the mismatch strain, we consider the layer 
as elastic with energy density 



E 



w(,e-eoj 



2(1 + J/) [l-2i^ 



-[tr(e-eo)]'+tr[(e-eo)2] 



We consider displacement fields u(a;, y) that are pe- 
riodic in X, u(l,y) = u(0,2/) + (const. )i for y G [0,ft.]; 
therefore, modulo a rescaling, we may restrict attention 
to a cell of unit length fi = [0,1] x [0,/i]. We assume 
that the layer cannot separate from the substrate, so that 
u(a;, 0) • j = for x G [0, 1] . Thus, letting u{x) = u(a;, 0) • i 
denote the tangential displacement of the layer at the in- 
terface, we define the incoherency strain ^{x) by 



l{x) 



du{x) 
dx 



and refer to the interface as coherent if 

7(a;)=0 a:e[0,l], 

and incoherent otherwise. 

We assume that the interfacial energy is an even, con- 
tinuous function /(7) of the incoherency strain, and that 
/(O) = 0, while /(7) > for 7 ^ 0. The total energy of 
the system is then given by 



J(u) 



w{e — Co) dxdy 



fh)dx. 



(1.1) 



II. CONVEXITY OF THE INTERFACIAL 

ENERGY DENSITY AND EXISTENCE OF 

SMOOTH EQUILIBRIUM STATES 

Equilibrium configurations of the layer correspond to 
states which minimize the total energy: classically, such 
configurations correspond to displacement fields u(x,y) 
such that 

J(u) = inin J(u) 
over a suitable^] space W of displacement fields on ft. 



with E and v Young's modulus and Poisson ratio, such 
that E > and -1 < 1/ < i 



^Assuming that / satisfies the growth condition f{-y) < 
C(l -I- I7I'), with C > and <? > 1, we minimize J on the 
space W of all functions u in the Sobolev space W^''^{Q,, TR^), 
whose restriction to dQ belong to W^'''{dQ, IR^), which satisfy 
the periodicity condition and are tangential to the interface. 



The chief problem here is that when / is not convex 
such mininiizers u may not exist. Even so, one can ob- 
tain valuable physical insight by studying minimizing se- 
quences; that is, sequences {u„} that tend to the infimum 
of the functional J(u) in the sense that 



^(u„) 



inf J(u), 
uew 



as n -^ cxD. In fact, as we shall see, such minimizing se- 
quences help to characterize the microstructures resulting 
from various choices of the interfacial energy /(7). 

Central to this method of analysis is the convex enve- 
lope /**(7) of /(7), which is the largest convex energy 
/**(7) with /**(7) < /(7) for aU 7; precisely, /** is the 
supremum of all convex functions g such that g < f. 
Then / is convex if and only if / = /** (Fig. 1). 




FIG. 1. The interfacial energy density (thick line) and its 
convex envelope (thin line). 

The convexity of /** allows us to obtain a unique so- 
lution u** of the "regularized problem" 



J(u**) 



Eh 



2(1-1/2) 



(7"-eo)' + /(7**), (2.3) 



then 



(«) /(7**) = r*(7**) ^ J(u**)= inf J(u), 

uGVv 

(b) /(7**)>r*(7**) ^ J(u")> inf J(u). 

uew 

Hence only when / is convex at 7**, does a minimum 
exist for the energy functional, at least in the classical 
sense. 



III. MINIMIZING SEQUENCES AND 
INTERFACIAL MICROSTRUCTURES 

Fix a thickness of the layer h > 0, let 7** be the average 
incoherency strain corresponding to the infimum of J, 
and assume that 7** ^ and /(7**) > f**{-f**), so 
that no smooth equilibrium state exists. In this case, 
physically meaningful results may be still be obtained by 
inspec tion of the minimizing sequences {u„}. 

By (2^) we may indeed restrict attention to the cor- 
responding sequences of incoherency strains 7„ at the 
interface, which must, in turn, be minimizing sequences 
of the interfacial energy (cf . |l| ) : 



J(u* 



inf 
uew 



wle — Cq) dx dy + / f**{'y)dx 
n Jo 



In fact, this solution has the explicit form 

**/ N ** . , ^i<^a - 7**) . 

u (a;,y :=7 a^i + ^7:; r-^yj, 

(l-i/) 



(2.1) 



and corresponds to the constant value 7** of the inco- 
herency strain 7 defined in Section IV below. 

Consider the case in which / is not convex at 7** . Here 
(cf. ||l|]) any minimizing sequence must satisfyjj 



u„ 



u** in n, and 



7„ dx -> 7* 



(2.2) 



as n — > 00. 

Thus minimizing sequences al ways converge in bulk 
to the homogeneous deformation {2A). As we shall see, 
minimizing sequences may not have a classical limit at 
the interface, but the "generalized limit"^, however it be 
visualized, corresponds to a well defined average inco- 
herency strain 7**. 

Of course, u** is a candidate minimizer of the total 
energy, but since 



hm / /(7„(a:))rfa; = r*(7**), 
/o 



(3.1) 



with J'**(^**) < /(7**). We shall consider two cases. 

• Oscillating sequences (Fig. 2). Consider the sit- 
uation in Figure 2(a), in which there is a A, < A < 1, 
such that 



7** = A7a -f (1 - A)7b, 

r*(7") = A/(7a) + (l-A)/(76). 



(3.2) 



j(-i) 














In 








lb 




~ 










~ 




7a 


- 








-! 









(a) 



(b) 



FIG. 2. (a) The interfacial energy density and its convex 
envelope; (b) a typical element of a minimizing sequence 7„. 



Then (p.l]) may be satisfied by a minimizing sequence 
7„ oscillating between 7a and jb on two subsets mixing 
finely as n ^ 00. More precisely, we may take 



■^Convergence of u„ is here in the sense of W^'^{fl, H'^). 
^ The generalized or weak limit of a sequence can be thought 



of as the limit of its "local integral averages" 



7" (a;) 



la 

lb 



fe k + XX 

p rfe+A fc+1 



), 



(3.3) 



with k — 0, . 
any n 



1 (cf. Figure 2(b)). Thus, since for 



/(7„) dx = A/(7a) + (1 - x)fhb) = rhn, 



then (3T) is satisfied. Indeed, the 7„ are minimizers of 
the intcrfacial energy functional for fixed 7**, but their 
weak hmit 7(2:) = 7** is not. 

Thus, in the hmit as n — > 00, the above sequence de- 
scribes an interfacial microstructure whose corresponding 
incoherency strain takes the values ja and 7^ on infinites- 
imal patches of length fractions A and 1 — A respectively. 



• Concentrating sequences (Fig. 3). Condition (3.2) 
may not be satisfied in the important case in which / is 
strictly concave (see also Fig. 6). We have 



f**i"f) = '^l?!? with m : 



/(7) 



7 — >+oo 7 



lim 



/**(t 




(a) 



(b) 



FIG. 3. (a) The interfacial energy density and its convex 
envelope; (b) a typical element of a minimizing sequence 7™. 

Fix /i > as above and assume that 7** > 0. Then 



(3.1) can only be satisfied by sequences 7„ which become 
unbounded on smaller and smaller sets, for instance 



Inix) 



7""n 




X e 



(3.4) 



for /c = 0, . 



1. In this case, (3T) is satisfied since 



f{ln{x)) dx = > m7 







Note that, since the 7„ tend to grow without bound on 
small sets, the corresponding displacements m„ at the 
interface tend to develop microscopic jumps; these may 
be identified with interfacial dislocations. 



IV. CRITICAL THICKNESS AND SMOOTHNESS 

OF THE CONVEXIFIED INTERFACIAL 

ENERGY DENSITY 

We now study the behavior of the minimizers when 
h varies. We show that, when /** is smooth at zero, 
the interface is incoherent (i.e., 7** ^ 0) for any /i > 0, 
while if /** is non-smooth at zero, there exists a critical 



thickness he such that for h < he the interface is coherent 
(so that 7** = 0). 

Cons ider first a smooth /**. Let ^(7) denote the right 
side of (2.3) when 7** is replaced by 7. Then j is smooth 
and convex and 7** is the unique solution of ^'(7) — 0. 
Since a direct calculation shows that j'(0) can only vanish 
when h = then, for h > 0, 7** ^ and the interface is 
incoherent. 

Assume now that /**, and thus j, is non-smooth at 
zero. Then the condition j'(0) — must be replaced 
by e [j!_(0), j+(0)], where j± are the left and right 
derivatives of j. This condition is in general satisfied by 
an interval of values for the thickness (cf. O), and thus 
the equilibrium interface is coherent for all h < he, with 
he a suitable critical thickness. 

In general, the convexified interfacial energy /** may 
be non-smooth at other values 7 of the incoherency 
strain: if this is the case, we assume that f**{j) = f{j), 
so that a regular solution exists to the minimization prob- 
lem. Proceeding as above, we see that the condition 
which assures that 7** = 7 is that G [j!_(7), j+(7)], 
and this condition again determines an interval for the 
thickness for which the incoherency strain remains fixed 
at the value 7(cc) = 7** = 7: the film remains 'glued' to 
the substrate. 



V. DISCUSSION 

We now turn to the analysis of specific forms of the 
interfacial energy density / and discuss the existence of 
smooth equilibrium configurations for the film, the for- 
mation of microstructures at the interface, and the crit- 
ical thickness for incoherency. Remark that the height 
of the film h is now allowed to vary, so that the average 
incoherency strain 7** also varies accordingly. 

• / convex and smooth (Fig. 4(a)). Since in this 
case 7(7) — /** (7), a regular solution exists and is given 
by u** in ( pl] ). In fact, the restriction of u** to the 
interface, namely u{x) = j**x, is a smooth minimizer of 
the interfacial energy functional. 




(a) (b) 

FIG. 4. Interfacial energy density: (a) smooth and convex; 
(b) non-smooth but convex. 

Thus the film is uniformly incoherently strained with 
respect to the substrate, but no fine structure appears. 
Moreover, since f**{'y) is smooth, the interface relaxes to 
incoherency for any thickness of the layer. 



• / convex but non-smooth at 7 = (Fig. 4(b)), 

Again f{-f) = /**(7), so that the homogeneous displace- 
ment u** in ( |2.l| ) is a minimizer of J, and no fine struc- 
ture develops at the interface. 

Now, since / is non-differentiable at 7 = 0, there exists 
a critical thickness for the transition to incoherency such 
that for h < he the interface is coherent, while for h > he 
the interface is uniformly strained with respect to the 
substrate. 

This form of interfacial energy might be appropriate to 
describe "glassy" interfaces between crystals with large 
mismatch. 

• / nonconvex and nonsmooth (Fig. 5). We assume 
here that / has pointed minima at values 7^ G {70 = 
0,±7i,±72, . . .} of the incoherency strain, so that the 
convexified interfacial energy /** is piecewise linear in the 
intervals (0,7i), (71,72),. . ., but non smooth at 7 == 7^. 



f(i) 



f"{i) 




FIG. 5. Nonsmooth scalloped energy / and its convex en- 
velope /**. 



7** = 0, a coherent interface. For values of h at which 
7** 7^ coherency is lost and fine microstructures appear 
at the interface. 

/**(7) 




FIG. 6. Nonsmooth strictly concave energy / and its con- 
vex envelope /**. 

More precisely, observe first that /** is non-smooth at 
zero, so that there exists a critical thickness he for inco- 
herency. W hen h > he, since 7** 7^ 0, the displacement 
field u** in (2.1) is no longer a solution; the minimizing 
sequences are concentrating, and energy is minimized by 
allowing the incoherency strain to become infinitely large 
on infinitesimal intervals: accordingly, the interfacial dis- 
placement Un tends to develop microscopic jumps, which 
we may interpret as interfacial dislocations. 

These energies might describe small-misfit epitaxial 
layers. 
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